Search for: All records

Machine learning methods are increasingly being employed as surrogate models in place of computationally expensive and slow numerical integrators for a bevy of applications in the natural sciences. However, while the laws of physics are relationships between scalars, vectors and tensors that hold regardless of the frame of reference or chosen coordinate system, surrogate machine learning models are not coordinate-free by default. We enforce coordinate freedom by using geometric convolutions in three model architectures: a ResNet, a Dilated ResNet and a UNet. In numerical experiments emulating two-dimensional compressible Navier–Stokes, we see better accuracy and improved stability compared with baseline surrogate models in almost all cases. The ease of enforcing coordinate freedom without making major changes to the model architecture provides an exciting recipe for any convolutional neural network-based method applied to an appropriate class of problems. 

Recent methods to simulate complex fluid dynamics problems have replaced computationally expensive and slow numerical integrators with surrogate models learned from data. However, while the laws of physics are relationships between scalars, vectors, and tensors that hold regardless of the frame of reference or chosen coordinate system, surrogate machine learning models are not coordinate-free by default. We enforce coordinate freedom by using geometric convolutions in three model architectures: a ResNet, a Dilated ResNet, and a UNet. In numerical experiments emulating 2D compressible Navier-Stokes, we see better accuracy and improved stability compared to baseline surrogate models in almost all cases. The ease of enforcing coordinate freedom without making major changes to the model architecture provides an exciting recipe for any CNN-based method applied on an appropriate class of problems. 

Abstract Many nucleosynthetic channels create the elements, but two-parameter models characterized byαand Fe nonetheless predict stellar abundances in the Galactic disk to accuracies of 0.02–0.05 dex for most measured elements, near the level of current abundance uncertainties. It is difficult to make individual measurements more precise than this to investigate lower-amplitude nucleosynthetic effects, but population studies of mean abundance patterns can reveal more subtle abundance differences. Here, we look at the detailed abundances for 67,315 stars from the Apache Point Observatory Galactic Evolution Experiment (or APOGEE) Data Release 17, but in abundance residuals away from a best-fit two-parameter, data-driven nucleosynthetic model. We find that these residuals show complex structures with respect to age, guiding radius, and vertical action that are not random and are also not strongly correlated with sources of systematic error such as $\log \left(g\right)$,T_eff, and radial velocity. The residual patterns, especially in Na, C+N, Mn, and Ce, trace kinematic structures in the Milky Way, such as the inner disk, thick disk, and flared outer disk. A principal component analysis suggests that most of the observed structure is low-dimensional and can be explained by a few eigenvectors. We find that some, but not all, of the effects in the low-αdisk can be explained by dilution with fresh gas, so that the abundance ratios resemble those of stars with higher metallicity. The patterns and maps we provide can be combined with accurate forward models of nucleosynthesis, star formation, and gas infall to provide a more detailed picture of star and element formation in different Milky Way components. 

Abstract The element abundance pattern found in Milky Way disk stars is close to two-dimensional, dominated by production from one prompt process and one delayed process. This simplicity is remarkable, since the elements are produced by a multitude of nucleosynthesis mechanisms operating in stars with a wide range of progenitor masses. We fit the abundances of 14 elements for 48,659 red-giant stars from APOGEE Data Release 17 using a flexible, data-drivenK-process model—dubbedKPM. In our fiducial model, withK= 2, each abundance in each star is described as the sum of a prompt and a delayed process contribution. We find thatKPMwithK= 2 is able to explain the abundances well, recover the observed abundance bimodality, and detect the bimodality over a greater range in metallicity than has previously been possible. We compare to prior work by Weinberg et al., finding thatKPMproduces similar results, but thatKPMbetter predicts stellar abundances, especially for the elements C+N and Mn and for stars at supersolar metallicities. The model fixes the relative contribution of the prompt and delayed processes to two elements to break degeneracies and improve interpretability; we find that some of the nucleosynthetic implications are dependent upon these detailed choices. We find that moving to four processes adds flexibility and improves the model’s ability to predict the stellar abundances, but does not qualitatively change the story. The results ofKPMwill help us to interpret and constrain the formation of the Galaxy disk, the relationship between abundances and ages, and the physics of nucleosynthesis. 

Abstract Discoveries of gaps in data have been important in astrophysics. For example, there are kinematic gaps opened by resonances in dynamical systems, or exoplanets of a certain radius that are empirically rare. A gap in a data set is a kind of anomaly, but in an unusual sense: instead of being a single outlier data point, situated far from other data points, it is a region of the space, or a set of points, that is anomalous compared to its surroundings. Gaps are both interesting and hard to find and characterize, especially when they have nontrivial shapes. We present in this paper a statistic that can be used to estimate the (local) “gappiness” of a point in the data space. It uses the gradient and Hessian of the density estimate (and thus requires a twice-differentiable density estimator). This statistic can be computed at (almost) any point in the space and does not rely on optimization; it allows us to highlight underdense regions of any dimensionality and shape in a general and efficient way. We illustrate our method on the velocity distribution of nearby stars in the Milky Way disk plane, which exhibits gaps that could originate from different processes. Identifying and characterizing those gaps could help determine their origins. We provide in an appendix implementation notes and additional considerations for finding underdensities in data, using critical points and the properties of the Hessian of the density. 7 7 A Python implementation of t methods presented here is available at https://github.com/contardog/FindTheGap . 

Abstract In the coming decade, thousands of stellar streams will be observed in the halos of external galaxies. What fundamental discoveries will we make about dark matter from these streams? As a first attempt to look at these questions, we model Magellan/Megacam imaging of the Centaurus A (Cen A) disrupting dwarf companion Dwarf 3 (Dw3) and its associated stellar stream, to find out what can be learned about the Cen A dark matter halo. We develop a novel external galaxy stream-fitting technique and generate model stellar streams that reproduce the stream morphology visible in the imaging. We find that there are many viable stream models that fit the data well, with reasonable parameters, provided that Cen A has a halo mass larger than M 200 > 4.70 × 10 12 M ⊙ . There is a second stream in Cen A’s halo that is also reproduced within the context of this same dynamical model. However, stream morphology in the imaging alone does not uniquely determine the mass or mass distribution for the Cen A halo. In particular, the stream models with high likelihood show covariances between the inferred Cen A mass distribution, the inferred Dw3 progenitor mass, the Dw3 velocity, and the Dw3 line-of-sight position. We show that these degeneracies can be broken with radial-velocity measurements along the stream, and that a single radial velocity measurement puts a substantial lower limit on the halo mass. These results suggest that targeted radial-velocity measurements will be critical if we want to learn about dark matter from extragalactic stellar streams.